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Abstract. The oxygen (O/H) and N/O abundance ratios along 
the bar of 16 barred spiral starburst galaxies are determined 
using long-slit spectroscopy. The abundance gradients and the 
spatial distribution of the ionized gas along the bar are used to 
understand the role of bars in starburst galaxies. 

The oxygen abundance gradients are steeper than in nor- 
mal barred galaxies, with a mean of -0.15 dex/kpc, while the 
intersects are low. This excludes the possibiUty that starburst 
galaxies in this sample are chemically evolved galaxies rejuve- 
nated by the effect of a bar. The nitrogen-to-oxygen abundance 
gradients are flatter than the oxygen ones with a mean of -0.05 
dex/kpc. But N/O intersects are high, which rules out the pos- 
sibility that a large quantity of gas was recently funneled by a 
bar toward the center of a young galaxy. There is no correlation 
between the abundance gradients and the bar properties, which 
suggests that bars did not influence the chemical evolution of 
these galaxies. Therefore, bars cannot be at the origin of the 
bursts in the nuclei of our sample galaxies. 

The oxygen and N/O abundance gradients are generally 
stronger in the bar than in the disk and are linked together by a 
linear relation consistent with a primary + secondary origin 
for the production of nitrogen. This can be fully explained in 
terms of star formation history in these galaxies. The gradients 
build up from the inside out, becoming stronger as the oxygen 
and N/O abundances increase in the bulge while staying low in 
the disk. This behavior is consistent with a simple Schmidt law 
relating the density of star formation to that of gas. 

In many of the sample galaxies, star formation occurs at one 
or both ends of the bar. The low level of chemical enrichment 
in these regions suggests that they recently experienced bar- 
triggered star formation: this is the only visible effect of bars. 
Our analysis shows that bars probably appeared very recently 
(a few 10^ years) in the starburst galaxies, which are relatively 
"young" galaxies still in the process of formation. 
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1. Introduction 

Massive galaxies play a fundamental role in the chemical evo- 
lution of the Universe, because they are the only systems 
which meet the physical conditions necessary for the build- 
up of heavy elements. Today, by using simple techniques and 
modest-size telescopes, it is possible to understand how the 
production of metals occurs in nearby massive galaxies and 
to relate this production to their formation and evolution pro- 
cesses. 

While still in its early stages, the study of chemical abun- 
dances in external galaxies has already produced important re- 
sults. One of the most interesting is the discovery that galax- 
i es in the Universe se em to follow a mass-metallicity relation 
( Zaritsky et al., 1994 ). Although we still have to understand the 



physical causes behind this behavior, there is no doubt that this 
phenomenon is a key parameter for understanding galaxy evo- 
lution. With this prospect in mind, we recently verified that the 
massive starburst nucleus galaxies (SBNGs) follow the same 
mass-metallicity relation as normal ones (Coziol et al., 1997a). 
Our observations not only confirm the universality of the mass- 
metallicity relation for a wide range of galaxies, but also sug- 
gest that its origin may be linked to the different paths followed 
by early and late-type spiral galaxies to form their bulge and 



disk (Coziol et al., 1998) 
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* Based on observations obtained at the 193cm telescope of Obser- 
vatoire de Haute-Provence, operated by INSU (CNRS) 



Another important phenomenon, unveiled by chemical 
abundance studies, is that the metallicity in the ionized inter- 
stellar medium of spiral galaxies decreases outward (e.g. Vila- 
Costas & Edmunds, 1992, and references therein). Various hy- 
potheses have been proposed to explain this gradient (Pagel, 
1989; potz & Koppen, 1992| ), but their verification is difficult. 
The reason is that many parameters and processes involved in 
the formation of an abundance gradient are still not well un- 
derstood, such as the initial mass function, the yield, the mech- 
anisms and time scales for the formation of halos, bulges and 
disks, as well as the possible infall or outfall of matter. Dif- 
ferent numerical simulations have shown, however, that one 
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important ingredient for producing abundance gradients is a 
non linear dependence of the star formation rate on gas density 



Table 1. Catalogue elements of the sample galaxies and ex- 
cerpts from our observing logbook 



( Wise & Silk, 1989; Gotz & Koppen, 1 992 ; Molla et al. , 1 996 ) 



It has also been suggested that once an abundance gradient is 
established, the presence of gas flows can either amplify or re- 
duce it ( [Edmunds, 1990[ potz & Koppen, 1992| ). 

Attempts to relate abundance gradients to other characteris- 
tics of galaxies, such as morphological type, luminosity or cir- 
cular velocity, have also been performed (Vila-Costas & Ed- 
munds, 1992; pye & Kennicutt, 1993|; [Zaritsky et al., 1994|), 



but no clear correlation has emerged. One interesting trend is 
however observed: barred spiral galaxies seem to have shal- 
lower abundance gradients than non-barred ones (Vila-Costas 



& Edmunds, 1992; Edmunds & Roy, 1993; Zaritsky et al.. 



1994; Martin & Roy, 1994). From a theoretical point of view, 
bars are exp ected to funnel g as from the outer parts to the nuclei 
of galaxies ( Noguchi, 1988). The fl ow along the bar can induce 
strong mixing ( Friedli et al., 1994 ) and reduce metallicity gra- 
dients. But star formation induced by the bar can also be strong 
enough to counteract dilution effects and increase the gradient 
in the inner parts of galaxies. Steep abundance gradients have 
indeed been reporte d in the inner parts of some barred galaxies, 
such as NGC 3359 ( fVlartin & Roy, 1995| ) and NGC 1365 (Roy 
& Walsh, 1998). 

In the present paper, we analyse the abundance gradients of 
oxygen (O/H) and nitrogen to oxygen ratio (N/0) in 16 barred 
starburst galaxies. Our goal is to establish if starburst galaxies 
follow the trends observed in normal galaxies. We also want to 
determine if there is a relation between the bar and the burst of 
star formation. We have recently proposed that SBNGs may be 
"young" galaxies still in their formation process (Coziol et al. 
1997a,1998). This possibility raises several new questions. Is 
the initial distribution of abundances in starburst galaxies flat, 
or have they already had enough time to establish steep gradi- 
ents? Is it possible to determine how the gradient builds up in a 
young galaxy and establish whether the presence of the bar is 
important? Could bars be younger in starburst than in normal 
galaxies? 

Our analysis innovates by studying the abundance gradients 
of two key elements, oxygen and nitrogen, which are assumed 
( [Coziol et al., 1999a ) to be produced by different types of stars 
and, thus, to be released into the interstellar medium on differ- 
ent time scales. This gives us a deeper insight into the processes 
of chemical evolution, because we can test the influence of star 
formation and of structures, like a bar, over a relatively long pe- 
riod of time (a few Gyrs). This advantage should in particular 
enable us to determine whether bars are young and if they play 
a major role in triggering or feeding the nuclear starbursts. 

2. The sample of barred starburst galaxies 

The 16 galaxies studied in this paper are a sub-sample of 144 
Markarian barred spiral galaxies studied spectroscopically by 



Contini ( |1996D and Contini et al. ( |1998[ ) (Paper III). This sub- 
sample is a selection of galaxies which are experiencing very 
intense bursts of star formation, as estimated by their high 
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Ha luminosity (L(Ha) > 10''° erg s ^) or large H/3 equiv- 
alent width (EW(H/3) > 30 A). Three Wolf-Rayet galaxies 
are included in our sample: Mrk 710, 712 and 799. Another 
Wolf-Rayet galaxy satisfying our selection criteria, NGC 6764, 
(Contini et al., 1997), was later added to the list; but, as we will 



show later, it turns out to be a LINER and cannot be included 
in our analysis. 

Some catalog elements for the sample galaxies and 
excerpts from our observing logbook are presented 
i n Table p]. The morpholog ies were taken in NED 



: 75 km s ^ Mpc ^), and the radius at 25 mag 
Both values were found in LEDA. The last two 



(http://nedwww.ipac.caltech.edu), or in LEDA (http://www- 
obs.univ-lyonl.fr) for Mrk 13 and 306. Table |l] also lists the 
distance D, estimated from the recession velocity (assuming 
H 

arcsec 

columns indicate the exposure times and position angles of the 
slit which were used during our spectroscopic observations. 

3. Spectroscopic observations and reductions 

Most of the spectroscopic observations were obtained in Octo- 
ber 1995 and June 1996 at the 1.93 meter telescope of Obser- 
vatoire de Haute-Provence. The CCD was a thinned 512x512 
Tektronix (pixel size 27/im). We used the Carelec spectrograph 



(Lemaitre et al., 199C ) with a spectral dispersion of 260 A/mm, 
which covers the spectral range 3600 to 7200A, at a resolution 
of - 7A. The slit was aligned along the bar During the first run, 
we took 90-minute spectra of most galaxies, with a slit width 
of 3.0". For flux calibration, we observed the standard stars 
G191B2B, Hih 600, HD 217086 and BD -h28° 4211 (M assey 
et al., 1988). In the second run, we took three 25-minute spec- 
tra of Mrk 799 and two of NGC 6764, with a sUt width of 2.8". 
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One standard star, HD 192281, was used for flux calibration 



(Masseyetal., 1988). The galaxies Mrk 710 and 712 were ob- 



served during previous runs at the same telescope, using similar 



instrumental settings (Contini et al., 1998). The average seeing 



during the nights was of the order of 1.2". 

All spectra were reduced with MIDAS, applying standard 
procedures: offset and flatfield corrections, sky subtraction, 
wavelength and flux calibration, airmass and galactic extinc- 
tion corrections and elimination of cosmic impacts. In order 
to study abundance variations along the bar, series of one- 
dimensional spectra were extracted from the two-dimensional 
spectra of each galaxy by averaging three contiguous rows 
along the slit. Observed emission-line fluxes and equivalent 
widths were then measured in each elementary spectrum. For 
the rest of the reduction we followed the procedure outlined in 



o 
o 

in 



o 



Contini et al. (1995) (Paper 1): the Balmer (Ha and H/?) absorp- 



tion contamination due to underlying stellar populations was 
corrected by adding 2 A to the equivalent width of the Balmer 
emission lines; the principal Balmer decrement (Ha/H/3) was 
used to estimate the reddening coefficient ch/3 and to calculate 
absolute dereddened fluxes. 

Our analysis may differ from other studies found in the lit- 
erature, because we measured the spatial variations of physical 
parameters only along one direction (along the bar), and not 
across the whole two-dimensional disk. This difference should 
be kept in mind when comparing our results with those of other 
studies. 

4. Properties of emission-line regions 




-1.0 0.0 
log([NII]X6584/Ha ) 

Fig. 1. Diagnostic diagram for the emission-line regions of the 
17 sample galaxies. The curved solid line separates Hll region- 
like objects {left) from active galaxies (right). The horizontal 
solid line separates Seyfert 2 nuclei (top) and LlNERs (bot- 
tom). The dotted curve denotes the theoretical model of disk 
Hll region of McCall et al. (1985). Most of the emission-line 
regions are of Hll type, except those located in the central part 
of NGC 6764 which are of LINER type (asterisks) 



A standard diagnostic diagram ( Baldwin et al 
& 



Osterbrock, 1987) was used to determine the dominant 
source of excitation of the emission-line regions in our sam- 
ple. Most of these regions fall in the domain of Hll region- 
like spectra (Fig. |l]), referred to hereafter as Hll regions. The 
only exceptions are two regions in Mrk 545, which have weak 
signal-to-noise ratios, and the inner region of NGC 6764. 
There has been some confusion in the literature about the na- 
ture of the activity in this last galaxy (see Contini et al. 1996 
for the various classifications). Our new spectra clearly show 
that the center of NGC 6764, as well as the circum-nuclear re- 
gions, are the seat of LINER activity. Because LINER is a non- 
thermal phenomenon (implying either an AGN or shocks), the 
metallicity in these regions cannot be determined by the meth- 
ods used in this paper. Consequently, NGC 6764 (and the two 
regions in Mrk 545) have been excluded from our abundance 
gradient analysis. 

Our sample of Hll regions spans a large range of excita- 
tion levels, from high-excitation spectra (log ([Olll]/H/?) > 
0.4), characteristic of Hll galaxies, to low-excitation spectra 



198 1 ; Veilleux 5. Determination of oxygen and nitrogen abundances 

5.7. Determination of oxygen abundances: general method 



The oxygen abundances are determined using empirical rela- 
tions between log(0/H) and the emission-line ratios and 
i?23, which are defined as: 



i?3 = ([Oiii]A4959 + [Oiii]A5007)/H/3 



(1) 



(log ([Olll]/H/5) < 0.4) typical of SBNGs (Coziol, 1996). The 



low-excitation Hll regions have [N II]A6584/Hq; ratios which 
are relatively high. This is a common property of SBNGs (Co- 



ziol et al., 1997b; Contini et al., 1998), which possess a slight 



overabundance of nitro gen as compared to Hll regions with 
comparable metallicity ( Coziol et al., 1999a ). 



i?23 = ([Oii]A3727 + [Oiii]A4959 + [Oiii]A5007)/H/3 (2) 

For our analysis we use the empirical relations first estab- 
lished by Edmunds & Pagel (1984) and quantified by Vacca & 
Conti (1992). In general, the above two relations yield simi- 
lar results. In some cases, however, the relation based on i?3 
gives higher oxygen abundances than the one based on i?23, 
because the extinction correction affects differently the two re- 
lations; a high extinction correction tends to decrease the ratio 
i?3, while increasing the ratio i?23. The consequence of this ef- 
fect is that when dust extinction is high, one overestimates the 
oxygen abundance using R-^ and underestimates it using i?23. 
For this reason, we will use both relations and adopt the mean 
value. Exceptions are Mrk 7 10 and 7 12, where [O ll]A3727 was 
not observed and only i?3 was used. 
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5.2. Determination of extinction in the center of galaxies 

In the center of half of the sample galaxies (Mrk 12, 307, 326, 
332, 545, 799, 898 and 1088), Balmer absorption features due 
to an underlying intermediate-age stellar population severely 
affect the measurement of the Balmer emission lines. This pre- 
vents us from determining the extinction coefficient ch/j which 
is needed for calculating the abundances. The standard correc- 
tion of 2A on the equivalent width is not adequate for these 
cases (because we still cannot measure H/3 correctly) and it 
is thus necessary to find a more suitable method to solve this 
problem. 

It has been known for some time that there is an empiri- 
cal relatio n between \og{Rz) and lo g([N ll]/[0 III]) (Alloin et 
al., 1979; Edmunds & Pagel, 1984). The advantage of using 
the [N II] and [O III] emission lines for determining the oxy- 
gen abundance is that they are not affected by underlying stel- 
lar populations. Such a method is ideal for analyzing the cen- 
tral regions of galaxies. However, we still have to correct these 
lines for dust extinction. As we will now show, this can be done 
empirically. 

According to the definition of the extinction coefficient, 
the observed and dereddened fluxes of the [Nll]A6584 and 
[Olll]A5007 lines (hereafter [Nil] and [Olll]) are linked by 
the relation: 



log[Nii]/[Oin] 



dred 



log[Nii]/[Oiii]„fc, -0.3XCH/3 (3) 



If one plots log([Nll]/[Oni])dred versus log([Nll]/[Olll])ohs 
for a sample of spectra, statistically large enough to include 
evenly distributed values of ch/3 between 0.0 and ^ 1.5, one 
expects the data to cover uniformly the region between the line 
of zero extinction and that of ch/3 ^ 1-5. We notice instead that 
the data gather along an oblique line of slope smaller than unity, 
which intersects the lines of constant ch/3 (see Fig. This is 
observed in two samples of Hll regions in normal galaxies (Mc- 
Call et al., 1985; IVila-Costas & Edmunds, 19931), as wefl as in 



our own sample of starburst galaxies (excluding the spectra af- 
fected by Balmer absorption). This correlation translates the 
fact that, on average, the larger the value of log([Nll]/[Olll]) 
the larger that of ch/s- The reason behind this behavior is that 
log([Nll]/[Olll]) is correlated with oxygen abundance (see 
Sect. 5J and Fig. ||), which in turn is correlated with the amount 
of dust and consequently with internal extinction (Heckman et 
al., 1998). 

This observed correlation provides us with an empirical 
way of estimating ch/3 from the observed fluxes of [N ll] and 
[Olll]. Using the sample of galaxies of McCall et al. (1985) 
and our own sample (excluding the spectra with Balmer ab- 
sorption), we find the following relation between the observed 
and dereddened line ratios: 



log([Nii]/[Oiii])<j,,^ = 0.931og([Nii]/[Oiii])„,, - 0.16 (4) 

The correlation coefficient of this relation is 0.99 and the un- 
certainty on log([Nll]/[Olll])rf,.ed is less than 0.1 dex. This re- 




-2.0 -1.0 0.0 1.0 

log([NII]6584/[0III]5007)„bs 

Fig. 2. Dereddened versus observed flux ratios of [N ll] A6584 
and [O III]A5007. The thin line is a 1 to 1 relation. A set of dot- 
ted parallel lines represents 6 ranges of ch/3 varying from 0.25 
to 1.50 (see text). The sample of normal galaxies of McCall et 
al. (1985) is used as an example {open circles). The thick line 
is the regression on the data 



lation win be used in Sect. 5.3 to calculate ([N ll]/[0 III]) 



)dred^ 



and thus i?3, in the center of galaxies affected by Balmer ab- 
sorption. 

Using the same samples, we find the following relation be- 
tween CH/3 and ([Nii]/[Oiii])o6s: 



CH/3 = 0.451og([Nii]/[Oiii])„fc, + 0.51 



(5) 



The correlation coefficient is about 0.95 and the uncertainty is 
about 0. 1 in ch/3 . This relation will be used in Sect. 5 .4 to dered- 
den the [O III] A5007, [O ll] A3727 and Ha fines in the center of 
galaxies with strong Balmer absorption, and hence to calculate 

-^23- 

5.3. Calculating R3 in the center of galaxies 

Foll owing an idea developed by various authors (A lloin et al., 
1979; |3dmunds & Pagel, 1984| ), we now search for an em- 
pirical relation between i?3 and [N ll]/[0 III] in our sample of 
16 barred starburst galaxies. For this analysis we use only the 
dereddened fluxes of Hll regions, about 300 of them, which 
are not affected by Balmer absorption. The two quantities ap- 
pear to be strongly correlated, as shown on Fig. and we fit a 
two-degree polynomial to the data, which yields the following 
relation: 



log(i?3 



0.123 (log([Nll]/[Olll]))2 
0.661 log([Nii]/[0 III]) + 0.086 



(6) 



with a dispersion in log(i?3) equal to 0.07 dex. 

To check the validity and the possible "universality" of the 
above relation, we now apply the same procedure to two sam- 
ples of Hll regions in normal galaxies: those of McCall et al. 
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-2.0 -1.0 0.0 1.0 

log([NII]6584/[0III]5007)d,gd 



Fig. 3. Relations between log(i?3) and dereddened 
log([Nll]A6584/[Olll]A5007). a for the Hll regions in 
our sample of starburst galaxies that are not affected by Balmer 
absorption, b for 278 Hll regions in normal galaxies, from the 
samples of McCall et al. (1985) and Vila-Costas et al. (1993). 
The two continuous lines in a and b are adjusted polynomials 




8.5 9.0 
12 + log(0/H) 



Fig. 4. N/O abundance ratio as a function of oxygen abundance 
for the Hll regions in the sample galaxies. The different lines 
indicate the relation for primary, secondary and primary + 
secondary (dashed curve) origin of nitrogen (Vila-Costas & 
Edmunds 1993). The solar value is indicated by a large double 
circle. The two dotted curves are envelopes corresponding to 
the observational errors. 

Figure ||a shows that our data do not reach oxygen abun- 
dances as high as in normal galaxies. This is precisely the do- 
main where we cannot determine in the standard way, be- 
cause, in starburst galaxies, such high values for oxygen abun- 
dance are observed only in the nuclear regions, which are af- 
fected by Balmer absorption features. Fortunately, the parallel 
empirical relation for normal Hll regions (Fig. |^) continues 
smoothly into the high oxygen-abundance regime. We there- 
fore consider it legitimates to use Eq.^j to calculate i?3 in this 
range of oxygen abundances. 



(1985) and of Vila-Costas & Edmunds (1993). As shown in 
Fig. the behavior is the same. Fitting a two-degree polyno- 
mial we obtain: 



log(i?3) = - 0.133 (log([Nii]/[Oiii]))2 

- 0.709 log([Nii]/[0 III]) - 0.008 



(7) 



The dispersion in log(i?3) is 0.12 dex, that is slightly higher 
than for our sample. 

It is important to note that the two fitted curves in Fig. ^ 
differ by a horizontal shift of 0.14 dex. This shift is consistent 
with the overabundance of nitrogen in SBNGs with respect to 



normal Hll regions (Coziol et al., 1999a). Consequently, there 



is no unique relation between log(i?3) and log([Nll]/[Olll]), 
but specific relations for homogeneous samples of galaxies, 
such as normal ones or SBNGs. For a given value of the ra- 
tio [Nll]/[Olll], a starburst galaxy will have a lower oxygen 
abundance than a normal galaxy, or, for a given abundance, a 
starburst galaxy will have a higher [N ll]/[0 III] ratio than a 
normal galaxy. 



5.4. Calculating R23 and 0/H in the center of galaxies 

Once we know ch/3 in the center of the sample galaxies (using 
Eq. ||), it is relatively easy to estimate i?23. First, we correct 
[O III] A5007, [O II] A3727 and Ha for reddening using cr/j and 
estimate the dereddened H/3 flux using the theoretical Balmer 
decrement: H/Sdred = Ha(ired/2.85. Then, the dereddened ra- 
tio i?23 is computed assuming [Olll] = 1.35 x [Olll]A5007, 
and the oxygen abundance O/H is deduced using the relation 
ofVacca&Conti ( |1992| ). 

Comparing the abundances obtained using R3 and R23 in 
the center of galaxies, we find the same kind of sensitivity to 
dust extinction as that observed in the outer parts of galaxies. 
Consequently, we also adopt the average value of the oxygen 
abundances estimated by i?3 and R23 in the center of our sam- 
ple galaxies. 



5.5. N/O abundance ratios 

With the data at hand, it is possible to derive the nitrogen-to- 
oxygen abundance ratio (N/O). For Hll regions in our sample. 
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Table 2. Observed gradients and intersects in O/H and N/O 



Galaxy O/H N/O 



Mi-k# 
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(dex/kpc) 


Intersect 
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-0.097 ± 0.025 






545 


-0 


014 ± 0.008 


+0.04 ± 


12 


-0 


248 ± 0.045 




0.000 




0.000 


-0.91 ± 0.08 


-0.035 ± 0.015 




0.000 


602 


-0 


113 ± 0.015 


+0.21 ± 


14 


-0 


114 ± 0.052 




0.000 


-0 


095 ± 0.016 


-0.68 ± 0.16 


-0.147 ± 0.033 


-0 


129 ± 0.048 


710 


-0 


119 ± 0.049 


+0.23 ± 


08 


-0 


119 ± 0.049 


















712 




0.000 


-0.40 ± 


06 




0.000 


















799 


-0 


086 ± 0.010 


+0.54 ± 


12 


-0 


116 ± 0.021 


-0 


062 ± 0.035 


-0 


039 ± 0.008 


-0.64 ± 0.09 


-0.036 ± 0.019 


-0 


098 ± 0.031 


898 


-0 


184 ± 0.010 


+0.55 ± 


02 


-0 


184 ± 0.010 








0.000 


-0.58 ± 0.14 


0.000 






1076 


-0 


131 ± 0.010 


+0.19 ± 


02 


-0 


131 ± 0.010 






-0 


038 ± 0.020 


-0.85 ± 0.05 


-0.038 ± 0.020 






1088 


-0 


061 ± 0.009 


+0.47 ± 


09 




0.000 






-0 


027 ± 0.006 


-0.65 ± 0.06 


0.000 







the N/O abundance ratio has been derived from the dereddened 
emission Unes [Oll]A3727, [Olll]A5007 and [Nll]A6584 fol- 



lowing the method of Thurston et al.(1996), as described in 
Coziol et al. (|1999a|). 



The N/O abundance ratio as a function of oxygen abun- 
dance is shown in Fig. ^. Also shown in this figure are the ex- 
pected relations for a primary, a secondary and a primary + 
secondary origin for the production of nitrogen (Vila-Costas 
& Edmunds, 1993). The Hll regions in our sample seem to fol- 
low the primary + secondary relation. But the evolution of 
N/O with oxygen abundance does not trace a continuous be- 
havior: it suddenly rises and forms a sort of plateau between 
8.7 < 12 + log(0/H) < 9.1. This behavior has been interpreted 
as evidence for chemical evolution by a sequence of bursts of 
star formation (Coziol et al. 1999a). 

6. Results 

6.1. Spatial distribution of oxygen abundances and of Ha 

The oxygen abundance and Ha emission as a function of po- 
sition along the slit are presented in Fig. ^. The emission peak 
is usually centered on the nucleus of galaxies. Exceptions are 
Mrk 12, 13, 712 and 799. In the case of Mrk 712 and 799, the 
maximum intensities correspond to very young star-forming 
regions containing numerous Wolf-Ray et stars. When studying 
the gas distribution inside the bar, we find two possible situa- 
tions: more than half of galaxies have ionized gas covering the 
whole bar, while in the rest of galaxies the ionized gas is mostly 
confined to the center. A large number of galaxies (62%) also 
show ionized gas either at the two ends of the bar (25%), or, 
more frequently, only at one end (37%). 

Examining the spatial distribution of the oxygen abun- 
dance, we find that it usually reaches a maximum in the center 
of galaxies. In general, this maximum corresponds to that of 
the Ha emission. In some galaxies, there is also strong nebu- 
lar emission outside the nucleus, accompanied by a local rise in 



oxygen abundance (e.g. Mrk 306, 332 and 545). In other galax- 
ies this secondary peak of emission is not accompanied by an 
increase in oxygen abundance (e.g. Mrk 12, 13 and 799), and 
in some galaxies the oxygen abundance stays high even with- 
out any secondary peak of Ha emission (e.g. Mrk 307, 373 and 
712). 

6.2. O/H and N/O abundance gradients 

To study the abundance gradients, we have folded the profiles 
of O/H and N/O abundance ratios about a center of symmetry. 
This center is the locus of the principal peak of Ha emission 
(see Fig. ||). 

To make the analysis easier, we adjusted lines by least- 
square fits to the data in two different regions of each profile: 
one includes the bar and the nucleus and another the region 
outside the bar. The latter region corresponds to the disk. Note, 
however, that the profiles never cover the entire length of the 
disk (that is out to a surface brightness of 25 mag arcsec"^), 
simply because we did not detect ionized gas that far. A global 
gradient is also estimated by fitting a line to the data over the 
full extent of the ionized gas distribution. The intersect (central) 
abundance is determined by extrapolating this global gradient 
up to the center. The two regressions are plotted as continuous 
lines and the global gradient as a dotted line in Figs. ^ and ^ 
The profiles on the two sides of the center are distinguished 
by different symbols. The numerical values of different gradi- 
ents and the intersects are given in Table ^. The uncertainties 
reported in that Table are those obtained from the regression 
calculations. 

The oxygen abundance profiles as a function of radius (in 
arcsec) are shown in Fig. ^. They show complicated patterns. 
Although the uncertainties on the oxygen abundance (^ 0.2 
dex) prevent us from interpreting the multiple inflections of 
these profiles, we believe that part of these features must be 
real, reflecting variations in oxygen abundance produced by 
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Fig. 5. The oxygen abundance and Ha emission of the galaxies versus position along the slit in arcsec. The oxygen abundance 
is given in solar units, adopting log(0/H) ~ —3.16 for the solar value (Grevesse & Noels 1993). For each galaxy, a vertical 
continuous line marks its geometrical center and two vertical dotted lines indicate the two ends of the bar. In Mrk 710 the ends of 
the bar fall outside the figure (at 33.7 arcsec). When the spatial Ha profile is not a pure gaussian, we also draw its decomposition 
into multi-gaussian components 



different intensities of star formation along the bar. One good 
example is Mrk 332, where there is an increase of oxygen abun- 
dance linked to star formation at the two ends of the bar. 



The abundance gradients (continuous and dotted lines in 
Fig. ^) also show various patterns. In Mrk 2, 373, 710, 712, 
898 and 1076, the gradient can be estimated in the bar only. In 
these cases, the bar gradient is considered the global gradient. 
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Fig. 6. Radial profiles of oxygen abundance log(0/H) in solar units. The radial coordinate is in arcsec, with the origin at the 
center of symmetry of the profile. Values on both sides of the center are distinguished by two different symbols. In each panel 
the vertical arrow indicates the end of the bar. The slopes of the lines adjusted by least-square fits are given in Table ^. The two 
continuous lines correspond to the bar and disk abundance gradients and the dotted line corresponds to the global gradient 
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Fig. 7. Radial profiles of nitrogen-to-oxygen abundance ratio log(N/0). The solar value on this scale corresponds to —0.9. The 
meaning of the symbols is the same as in Fig. S The slopes of the lines adjusted by least-square fits are given in Table 
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Half of these gradients are relatively shallow —0.03) and 
half are moderately strong —0.15). 

For galaxies with ionized gas also in the disk, we distin- 
guish two cases: the disk abundance gradient can be either flat 
or negative. In Mrk 306, 326, 545 and 602, it is flat. Consider- 
ing that norma l spiral galaxies usua lly have monotonous nega- 
tive gradients ( Zaritsky et al., 1994 ), this makes these galaxies 
with flat disk gradient somewhat peculiar, although shallower 
outer profiles have been ob served before, for instances in NGC 
1365 dRoy & Walsh, 1998|), NGC 33 19 dZaritsky et al., 1994| ) 
and NGC 3359 ( [Martin & Roy, 1995| ). Other galaxies with ion- 
ized gas in the disk look more "normal", with negative disk 
gradients. In general, the gradients are shallower in the disk 
than in the bar There are two exceptions, Mrk 13 and Mrk 332. 
In the case of Mrk 332, the intense star formation located at the 
ends of the bar is obviously responsible for the steeper gradi- 
ent in the disk (see Fig. ||). For Mrk 13, there is no gradient in 
the bar simply because this bar is very short, it is the shortest 
bar among the 125 ones measured by Chapelon et al. (1999) 



(paper V). The bar of Mrk 1088 also seems to have a flat gradi- 
ent (Fig. ||), but the unfolded oxygen abundance profile (Fig. ^ 
shows that there is in fact a steep gradient accross the center 

In general, our global abundance gradients are negative. 
This reflects the fact that the oxygen abundance is almost al- 
ways highest in the center of the galaxy. One exception is Mrk 
712, which does not have any gradient at all. 

We determined the N/O abundance ratios of all galaxies ex- 
cept Mrk 710 and 712 (see Sect. The N/O profiles are 
shown in Fig. 0. In general, they look similar to those of oxy- 
gen. The N/O abundance ratio also increases toward the nu- 
cleus of the starburst galaxies, but the gradients are signifi- 
cantly shallower and many more galaxies have zero gradient. 
Further comparisons of the two abundance gradients are post- 
poned to Sect. |7| Values for the different gradients and the in- 
tersects are given in Table 0. 



7. Analysis of the abundance gradients 

7.1. Oxygen abundance gradients 

The presence of a bar could affect the abundance gradients in 
our sample of galaxies in two ways. By funneling metal-poor 
gas from the outer regions toward the center, the bar is expected 
to dilute the chemical elements, and thus reduce any initially 
present radial abundance gradient. But a bar could also stimu- 
late star formation at the centre, which, if strong enough, could 
establish or maintain a steep abundance gradient. In the follow- 
ing analysis, we look for such possible effects of the bar in our 
sample of starburst galaxies. 

Of the three categories of oxygen abundance gradients esti- 
mated here, the global gradient is the one which is most compa- 
rable to those from other studies. We show in Fig. ^the relation 
of global gradients and intersects with absolute magnitudes. 
This Figure should be compared with Fig. 1 of Edmunds & Roy 
(1993) for a sample of normal galaxies. The barred starburst 
galaxies generally have steeper gradients than normal barred 
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Fig. 8. a Global oxygen abundance gradient as a function of 
absolute magnitude. Galaxies where the global gradient is also 
the bar gradient are indicated by open squares; filled circles 
represent the other cases. The dashed line marks the mean gra- 
dient of normal unbarred galaxies with absolute B magnitude 
between-19 and -22 in the sample of Edmunds & Roy (1993). 
The continuous curve is a linear regression on the filled circles 
(correlation coefficient — 86%). b Intersect as a function of the 
absolute magnitude. The meaning of the symbols is the same 
as in a. The dashed line indicates a mean value of the central 
oxygen abundance in normal unbarred galaxies with absolute B 
magnitude between -19 and -22 (from Edmunds & Roy 1993) 



galaxies (Edmunds & Roy, 1993; Zaritsky et al., 1994). On av- 
erage, the gradients in the barred starburst galaxies are even 
steeper than those measured by Edmunds & Roy (1993) in nor- 
mal unbarred galaxies of comparable luminosity. This shows 
that the oxygen abundance in the sample of barred starburst 
galaxies is not diluted by the mixing effect of the bar. 

In the above analysis, six of the global abundance gradi- 
ents are in fact bar gradients, because no ionized gas was ob- 
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served in the disks of these galaxies. We should therefore re- 
examine Fig. 1^ and distinguish these cases from those where 
the global gradient includes the disk gradient. With this dis- 
tinction, we find that the global gradient in Fig. ^ tends to 
steepen as the absolute magnitude increases. A linear fit per- 
formed on the data yields a correlation coefficient of 86%. But 
if we exclude the extreme point (Mrk 13) at Mb ~ —18, this 
coefficient drops to 56%. This is more a trend than a true corre- 
lation. Galaxies for which the global gradient is the bar gradient 
also show a similar trend, but slightly displaced towards lower 



values. An explanation for this trend is proposed in Sect. |8^ 
Examining now the central oxygen abundances (Fig. 



we find that they are on average lower in barred starburst galax- 
ies than in normal unbarred galaxies. Because the global abun- 
dance gradients in barred starburst galaxies are strong, these 
low central abundances cannot be attributed to an artifact of the 
bar, that would have lowered the central abundance of chemi- 
cally evolved galaxies by dilution. This means that the starburst 
galaxies in this sample have lower oxygen abundances than 
normal ones. This is in agreement with our previous finding 
dCoziol et al., 1997a| ; |Coziol et al., 1998| ): SBNGs are chemi- 
cally less evolved than normal galaxies with comparable lumi- 
nosity and morphology. We h ave shown elsewhere (Coziol et 
al., 1997a; Coziol et al., 1998), that the only way to explain this 



phenomenon is by assuming that SBNGs are "young" galaxies. 
7.2. N/0 abundance gradients 

We have shown that trends in the oxygen abundance gradients 
are inconsistent with the mixing effect of a bar, unless the bar 
also triggers a central starburst which more than compensates 
for the dilution. However, the fact that we find both strong oxy- 
gen abundance gradients and lower intersects than in normal 
galaxies rules out a bar-triggered burst in the center of an "old" 
galaxy. 

But the bar could be responsible for starbursts in the cen- 
ter of a relatively "young" galaxy. In this case, there are two 
alternatives. The gas funneled by the bar could have been fed 
directly into the central region in a short time (a few Myrs), 
prompting a new burst of star formation that is sufficiently 
strong to produce steep abundance gradients. The other alter- 
native is that the bar acts over a much longer time scale (a few 
Gyrs), progressively feeding the inner regions with gas that is 
immediately recycled (since no mixing effects are seen). 

We can test the above two alternatives by studying the 
N/O abundance gradient s and c omparing them with the oxy- 
gen ones. In Coziol et al.( 1999a ), we have shown that the bulk 
of nitrogen in SBNGs is probably produced by intermediate- 
mass stars, whereas oxygen is mainly released by massive 
stars. In this scenario, the time scale for enrichment in nitro- 
gen (400 Myr < t < 2 Gyr) is much longer than the one for 
oxygen (r < 20 Myr). This allows us to look for effects of the 
bar on a time scale covering a few Gyrs. 

In Fig l^a, we show the relation between the global N/O 
abundance gradient and absolute magnitude. In this Figure, we 
distinguish galaxies with a gradient measured only along the 
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Fig. 9. a Global N/O gradient as a function of absolute mag- 
nitude. The symbols have the same meaning as in Fig. ^. The 
continuous line is a regression on the N/O values, with a corre- 
lation coefficient of 75%. For the sake of comparison, we also 
give the linear regression (dotted line) obtained for the gradient 
in O/H (see Fig. b Intersect as a function of absolute magni- 
tude. The horizontal dashed line corresponds to the solar value 
of N/O 

bar from those where the gradient also includes the disk. We 
observe the same trend for N/O as for oxygen: galaxies with 
a lower luminosity seem to have a stronger negative gradient. 
The relation implied by this tendency has a weak correlation 
coefficient (60% with Mrk 13 excluded) and a shallower slope 
than that found for oxygen. Again, galaxies with only a bar 
gradient seem to have shallower gradients than the other ones. 

Despite the small gradients, the central N/O values, shown 
in Fig. are rather high compared to the solar abundance. 
This implies that many generations of intermediate-mass stars 
have already contributed to enrich these galaxies in nitrogen, 
or, equivalently, that enough time has passed in the chemical 
evolution of these galaxies. 
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Why is the gradient systematically weaker for N/O than for 
oxygen? If this is due to the bar, it means that the effect of 
the bar appears on a time scale comparable to the nitrogen en- 
richment phase. This is unlikely as it implies that the oxygen 
enrichment phase, and the ensuing strong oxygen abundance 
gradient, happened after the nitrogen one. Therefore, as a mix- 
ing mechanism, the bar should have affected both elements in 
the same way, which is not observed. Our explanation for the 



observed difference is given in Sect 7.3. 

The fact that the N/O intersects are high rules out the possi- 
bility that the bars recently triggered central starbursts in these 
galaxies. A recent burst would have increased the oxygen abun- 
dance and decreased the N/O abundance ratio. The stronger the 
burst, the stronger the decrease in N/O. We would thus expect 
galaxies with the steepest oxygen abundance gradients to have 
the flattest N/O ones. In Fig. |lO^, we show the relation between 
the oxygen and N/O abundance gradients. The behavior is con- 
trary to what is expected. The bars consequently did not re- 
cently trigger bursts in these galaxies. 

73. The progressive build up of abundance gradients 

What is then the reason for the difference between the two 
gradients? In fact, the N/O abundance gradient is not inde- 
pendent of the oxygen one. This is shown in Fig. [l0| a. A lin- 
ear regression on the data yields the relation: A(N/0) / Ar = 
0.53 (±0.10) A(0/H)/Ar + 0.00 (±0.01) with a correlation 
coefficient of 89%. Surprisingly, the slope of this relation is al- 
most the same as the one found between the N/O and oxygen 
abundances: a linear fit performed on the data in Fig. ^ yields 
log(N/0) = 0.551og(O/H) +0.8 (|Coziol etal., 1999a[), which is 



consistent with a mixture of primary + secondary origin for 



nitrogen ( McGaugh, 1991 ) 



The relation between the two gradients is thus explained 
by the behavior of the oxygen abundance in the nuclei of the 
galaxies. This is shown in Fig. where the oxygen abun- 
dance gradient steepens when its intersect increases. This rela- 
tion is trivial if we assume that all disks initially have about the 
same low oxygen abundance (see Fig. Since the nitrogen 
enrichment is related to the increase in oxygen (See Sect. |5.5| 
and Fig. Q), the N/O abundance gradient grows with that of 
oxygen. The origin of both gradients in starburst galaxies can 
thus be explained solely by their star formation histories. 

8. The role of bars in starburst galaxies 

The only way to reconcile the above observations of abundance 
gradients with the possible effect of a bar is by assuming that 
the latter is a slow process; over a few Gyr period, the bar pro- 
gressively feeds the nucleus with gas which is immediately re- 
cycled into stars. This scenario requires a very fine tuning be- 
tween the feeding of the gas in the center of galaxies and star 
formation. It also requires that bars are relatively stable over a 
long period of time (a few Gyrs). We can test if these condi- 
tions apply to galaxies in our sample by looking for relations 
between the bar properties and star formation. 
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Fig. 10. a Relation between global gradients in N/O and O/H. 
The symbols have the same meaning as in Fig. ^. The contin- 
uous curve is a linear regression with a correlation coefficient 
of 89%. b Relation between the global gradients in O/H and 
its intersect. The continuous curve is a linear regression with a 
correlation coefficient of 68% 



8.1. Bar properties and star formation: evidence for young 
bars 

Numerical simulations predict that a strong bar maintains a 
steeper abundance gradient in the bar than in the disk for about 
1 Gyr, then both become comparably low (~ —0.02), and that 
the bar and disk gradients are never very different in the pres- 
ence of a weak bar (|Martinet & Friedli, 1997). We show in 



Fig. ^ that almost all galaxies in our sample show a stronger 
gradient in the bar than in the disk. Only two galaxies, Mrk 13 
and 332, have significantly steeper gradients in the disk than 
in the bar, but this peculiarity can be explained by other phe- 



nomena (see Sect. 5.2). In other words, the difference between 
bar and disk gradients cannot be due to an old bar It cannot 
be due to the effect of a young bar either, because, as we pre- 
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viously noted, the N/O abundance ratios would be low and the 
gradients would be anti-correlated. 

A simple explanation for the stronger gradients in the bar 
than in the disk is that the volume density of star formation 
is proportional to the gas volume density to some power n 
(1 < n < 2), as proposed by Schmidt's law (1959, 1963). 
Assuming that the chemical enrichment is proportional to the 
star formation, a spherical density distribution yields a vari- 
ation with radius in i?^'^/". The gradient measured from the 
center to the end of the bar (or to any intermediate radius) will 
necessarily be steeper than the one measured further out. 

The presence of ionized gas is a tracer of star formation in 
galaxies. Half the galaxies in this sample show the presence of 
ionized gas over the whole bar, while in the other half, the gas 
is mostly concentrated in the inner parts of the bar and in the 
nucleus. If star formation cannot result from gas flow toward 
the inner regions by the dynamical action of a bar, could it be 
linked to any other property of the bar, such as its strength or its 
length? The bar strength is identified with the deprojected bar 
axis ratio (b/a)^, and its length is given by the ratio of the de- 
projected major axis to corrected blue isophotal diameter 2a/Dc 
dChapelon et al., 1999[ ) (Paper V). A bar with a ratio b/a < 0.5 
is considered "strong" and a bar with a ratio 2a/Dc > 0.2 is 
considered "long". We also identify separately galaxies with 
various fractions of ionized gas along the bar and galaxies with 
early- (earlier than SBbc) or late-type (SBbc or later) mor- 
phologies. 



As shown in Fig. 12, most galaxies in this sample have a 
strong bar, which is a general charac teristic of the sample o f 
Markarian barred starburst galaxies ( Chapelon et al., 1999 ). 
There is a trend for galaxies with a longer bar to have ionized 
gas over a smaller extent of the bar length, regardless of their 
morphology. This is not an observational selection effect, be- 
cause three galaxies do not follow the above trend in Fig. |l2[ 
namely Mrk 306, 712 and 1076. In the case of Mrk 306 and 
1076, a long and strong bar with ionized gas over its whole 
length may be explained by the fact that they are both inter- 
acting with close companions (the rest of our galaxies are iso- 
lated). Mrk 712 looks like an isolated galaxy, but with a pecu- 
liar morphology. The peak of star forming activity is signifi- 
cantly displaced from the center of the galaxy (see Fig. |]). This 
active region may be ass ociated with a second nucleus ( Maz- 
zarella & Boroson, 1993), which makes this galaxy a strong 
candidate for a recent merger. A simple relation with the length 
of the bar thus explains the radial distribution of ionized gas; 
star formation does not cover the whole bar just because it is 
too long. 

We have shown so far that bars, despite the fact that they are 
strong, cannot account for the observed abundance gradients. 
This may be surprising in view of the convincing evidence for 
dynamical effects of strong bars on star formation and abun- 
dance gradients in galaxies provided by numerical simulations. 
The only explanation for such a situation is that bars in the cur- 
rent sample of galaxies are too young (a few 10^ years) to have 
had any effect. 
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We verified that there are no relations between the abun- 
dance gradients and the bar properties. The oxygen and N/O 
abundance gradients in the bar and in the disk show the same 
behavior. No relation is found between the abundance gradi- 
ents and the bar strength or length. Nor is any relation found 
between the abundance gradients and the different concentra- 
tions of gas along the bar In general, therefore, the bars have 
not influenced the two gradients. This supports our assumption 
that the bars are young. 

8.2. Bar-induced star formation at the bar ends 

If we compare the distribution of ionized gas (i.e. star forma- 
tion) with the oxygen abundance, we find that the behaviors 
of these two parameters along the bars are frequently uncor- 
rected (see Fig. ^ and Sect. xA). Intense star formation does 
not necessarily enhance the oxygen abundance, as we can see, 
for example, in the extranuclear star forming regions of Mrk 
12 and 799. These regions must be too young (a few Myrs) 
and their massive stars did not had enough time to change the 
chemical abundance of their environment. When we find, on 
the other hand, a region filled with ionized gas, but where the 
abundance is significantly higher, we conclude that this must 
be a region where star formation was "stable" over a longer 
period of time. The galaxy nuclei in our sample are obvious lo- 
cations for such stable star-forming regions. The relatively high 
N/O abundance ratios and the high proportion of intermediate- 
mass stars mixed with ionized gas indicate that star formation 
persisted in these regions for quite a long time (a few Gyrs). 

We found no evidence that bars triggered the starbursts ob- 
served at the center of the sample galaxies. On the other hand, 
they may have induced star formation at the ends of the bar (e.g. 
Mrk 307, 332), and more frequently at only one end (e.g. Mrk 
12, 13, 306, 545 and 799). This phenomenon may be due to gas 
compression in these regions (Roberts et al., 1979). Star for- 



mation induced by the bar at its ends has increased the oxygen 
abundance there. The result is a shallower abundance gradient 
in the bar and a stronger one in the disk. This may explain why 
the global oxygen and N/O abundance gradients are weaker 
when they are measured only in the bar Chemical enrichment 
produced by star formation at these points being lower than 
what is observed in the nucleus, this suggests that fewer star 
formation episodes occurred there, or that they happened rela- 
tively recently. 

Among all the possible effects of the bar on the evolution 
of starburst galaxies, therefore, the only plausible one seems to 
be star formation induced at one or both ends of the bars. The 
distribution of star formation in this sample of starburst galax- 
ies is consistent with the assumption that bars are too young (a 
few 10^ years old) to be at the origin of the bursts in the galaxy 
nuclei. 

9. Summary and conclusion 

One important result of our analysis is that bars have not had a 
great impact on the star formation history and chemical evolu- 



tion of starburst galaxies. The most straightforward explanation 
is that bars are too young. Our observations also clearly show 
that these young bars cannot be at the origin of the nuclear star- 
bursts. The high N/O abundance ratios together with the high 
level of star formation in the nuclear region of galaxies suggest 
that star formation proceeds almost steadily (or as a sequence 
of bursts) over a few Gyr period. 

The fact that we find strong O/H gradients while the oxy- 
gen abundance is low in the galactic nuclei supports our inter- 
pretation that SBNGs are young galaxies stiU in their process 
of formation (Coziol et al. 1997a; 1997b). It seems that these 
starburst galaxies are still building their abundance gradients. 
This process can be fully explained in terms of their star for- 
mation history. The gradients build up from the inside out, be- 
coming stronger as the oxygen and N/O abundances increase in 
the bulge while staying low in the disk. This behavior is con- 
sistent with a simple Schmidt law relating the density of star 
formation to that of gas. 

According to these results, starburst galaxies seem to form 
their bulge first, and then their disk. If the disk is younger than 
the bulge and the bar forms in the disk, it is not surprizing that 
we find mostly young bars in starburst galaxies. The above sce- 
nario may also explain why our results are not consistent with 
predictions made by different models of bar formation. The 
initial conditions assumed in these models are far from those 
observed in young starburst galaxies. A bar forms in a galaxy 
which already has a large bulge but a disk which is probably 
still in formation. 

When star formation stops in the bulge and increases in the 
disk, as in normal spiral galaxies, the metallicity of the disk 
grows, decreasing the gradient. If the end product of a starburst 
galaxy is a normal spiral galaxy, the latter, therefore, should 
have lower metallicity gradients than the former. If they are 
stable, bars may also play a more important role as galaxies get 
older 

But how general is our conclusion? Is what we observe a 
trait of SBNGs? What is the relation between the bulge and disk 
formation in starburst galaxies and in normal spiral galaxies? 
These are the questions we address in our companion paper 
(Coziol et al. 2000). 
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